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We discuss the universal scaling laws of order parameter fluctuations in any system 
in which the second-order critical behavior can be identified. These scaling laws 
can be derived rigorously for equilibrium systems when combined with the finite- 
size scaling analysis. The relation between order parameter, criticality and scaling 
law of fluctuations has been established and the connexion between the scaling 
function and the critical exponents has been found. 



1 Introduction 

Averaged quantities are still commonly used for the description of many com- 
plex processes in physics, chemistry or astrophysics. This kind of approach 
is inprinted in our minds by an unquestionable success of thermodynamics 
in which fluctuations around the thermodynamic value are small in large sys- 
tems. It is our experience that the gross measures in, e.g., the particle reaction 
data are not sufficient to discriminate among models unless supplemented with 
more fine grained information, especially fluctuations and correlations of vari- 
ous kinds. Recent resurgence of interest in fluctuations in the strong interaction 
physics is due to the experimental possibility to measure event by event many 
strongly fluctuating quantities, e.g., the multiplicity of produced particles in 
ultrarelativistic collisions of leptons, hadrons, nuclei, or the charges (masses) 
of fragmentation products of highly excited heavy-ion residue. The immedi- 
ate question is then what is the information contained in fluctuations and, 
in particular, in the order parameter fluctuations , how do these fluctuations 
scale with the system-size and what is the relation between criticality and the 
probability distribution of order parameter fluctuations. 
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2 Order parameter fluctuations 

Let us suppose that the thermodynamic free energy of finite system depends 
on three parameters : r\ (the intensive order parameter) n e (the distance to 
the critical point) and N (the size of the system) . Widomld has proposed that 
close to the critical point, the free energy density in the thermodynamic limit 
scales as : 

/ (A^,Ae)^A 2 -«/ (r ? ,e) , (1) 

where a, (3 are usual critical exponents. There is no critical behavior in a finite 
system. However, a finite system behaves like an infinite one if the correlation 
length £ becomes comparable to the typical length L of the system. This 
is basically the argument of Fisher and Barber □ leading to the finite-size 
scaling analysis of critical systems. The pseudocritical point for a finite system 
appears at a distance e ~ cN^ 1 / 1 " 1 from the critical poinicl , where c is some 
dimensionless constant which can be either positive or negativcu. One can then 
deduce the scaling of critical free energy density at this point : 

f c (r,,N)^r, 2 -^4>(vN^) . (2) 

Assuming the hypcrscaling relation : 2 — a = vd , and using Rushbrooke 
relation between critical exponents : a + 2/3 + 7 = 2, one can write the total 
free energy F(ij, e, N) = Nf (rj, e) at the pseudocritical point as follows : 

F C ( V ,N)~ f (r)NTm,c) . (3) 

The canonicaLprobability density P[rj] to get some value of the order parameter 
r\ is given bjcl : 

P[ V ] = Zn- 1 exp(-prF(Ti, e, N)) , (4) 

where j3r is the inverse of temperature. Using Eq. ([3]), one obtains the 
partition function : 

Z N - N'^TW ~< |^| > . (5) 

It is then easy to see that the probability density P[rj] obeys the first scaling 
law : 

<\v\>P[ri\=*{zw) = *( r, ~ < : ] ?l > 

V < \v\ > 

= a(/3 T ) exp \ -0rf o ( - ^ - , c ) ) ( (i i 
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with the constant coefficient (3t independent of 77 , and the scaling function 
$(zm) depending on a single scaled variable : z^ = (77— < \r]\ >)/ < \rj\ > . 
The scaling limit is defined by the asymptotic behaviour of P[rj\ when 77 — > 00, 
< |ry| >^ 00, but (r)/< \rj\ >) has a finite value. The temperature-dependent 
factor a(/9r) is determined by the normalization of P[t]\ . One may notice that 
the logarithm of scaling function (^) corresponds to the non-critical free energy 
density at the renormalized distance e = c from the critical point. If the order 
parameter corresponds to the cluster multiplicity, like in the fragmentation - 
inactivation binary (FIR] procesflQ, then (|) can be written in an equivalent 
form to the KNO scaling?! , proposed some time ago as the ultimate symmetry 
of S - matrix in the relativistic field theorp . Relation between the KNO scaling 
and the phase transition in Feynman-Wilson-aas as well as the criticality of 
self-similar FIB process was studied as welMIll . Defining the anomalous 
dimension for an extensive quantity -ZV77 as : 

g = \hji N ^ OQ g N = \im N ^ 00 —^— (\n < N\i]\ >) , (7) 

a In N 

one can see that due to (||), the scaling @ holds when g = (7 + /3)/(7 + 2/3) . 
Consequently, g is contained between 1/2 and 1. Whenever the cluster-size 
can be reasonably defined for the second-order transition, like it is the case in 
percolation, Ising model or Fisher droplet model, the exponent r of the power- 
law cluster-size distribution : n(k) ~ kr T , satisfies additional relation£3 : 
7 + (3 = l/a and 7 + 2/3 = (r - l)/a , which yield : g = 1/(t - 1) . This 
means that r has to be contained between 2 and 3 when (^) holds. 

What may happen if the order parameter is not known exactly? To illus- 
trate this point, let us consider a quantity : m = N K — N-q , where 77 is the 
true order parameter and n is larger than the anomalous dimension g . For 
large N , |m| is of order N K . Writing (H) with m instead of 77, and taking into 
account : P[r)]dr] = P[m]dm , one finds the delta - scaling : 

<H ^PH^^hJ ^I^ ) , 8= 9 -<l , (8) 

\ < \m\>° J k 

with the scaling function ^(z^)) depending only on the scaled variable : 
z (8) = ( m — < \ m \ >)/ < \ m \ > S ■ According to (||) and (g), the logarithm of 
scaling function : 

ln$0(,5)) = -f3 T fo{z ( s),c) , (9) 

is directly related to the non-critical free energy , in either ordered (c > 0) 
or disordered (c < 0) phase. The relation 6 = g/n in (||) singularizes impor- 
tance of an extensive variable : m = N(l — 77) = Nf/ . 77 can be useful in 
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phenonicnological applications and plays an important role in the percolation 
studie£3 . One finds for this choice : < m >~ N , with algebraic finite-size 
corrections, and the delta-scaling (||) with 5 = g. Hence, P[Nfj\ allows to 
determine the anomalous dimension g and, consequently, the ratio of critical 
exponents (3 and 7 . 

Let us suppose now that the extensive parameter m is not critical , i.e., 
either the system is in a critical state but the parameter m is not critical, or 
the system is outside of critical region. The value of m at the equilibrium is 
obtained by minimizing the free energy. Let us suppose that the free energy 
F is analytical in the variable m close to its most probable value m* , i.e., : 

F ~ N^^(m — m*)* +1 . (10) 

In most cases 0=1, though, in general, 4> can take any positive integer 
value. Using ( p^| ) one obtains : < |m| >^ fi* 'N , where fi* is a positive (finite) 
number independent of N , and : 

4, 4, 
Z N ~ AU+i ~ < \m\ . (11) 

The probability density P[m] verifies the generalized scaling law (|J) : 

< \m\ > 5 P[m] = = exp [ f ^-<\m\>) \ \ ^ 

y V < \ m \ >* +1 / / 

but now 8 (= 4>/(4> + 1) < 1) is constrained by the value of <j> ..In the generic 
case <j) = 1 , 5 equals 1/2 and the scaling function is GaussiaiO. The second 
scaling ( |l2|) holds for < m >^ N but now with the exponential finite-size 
corrections. 

3 Results 

The above results apply to any second order transition, and , in particular, they 
are not limited to the Landau- Ginzburg theory of phase transitions. As an 
illustration, let us look at the bond percolation model ( Fig. 1). The intensive 
order parameter is the normalized mass of largest cluster : rj — S max /N . Fig. 
la (the upper left plot) shows the scaling function in 3D-percolation 

at the value of bond activation probability (p cr = 0.2482) corresponding to 
the critical activation in the infinite system. Results for different lattice sizes 
are superposing well in the scaled variables (^j) for 5 = 1. One should recall 
that is also a 'portray' of the free energy density at a renormalized 

distance from the critical point. The probability distribution P[Nrj\ changes 
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Figure 1: The order parameter fluctuations in the bond percolation model in 3D are 
calculated in the lattices : N = 10 3 and N = 14 3 , and plotted in the scaling variables of 
delta-scaling (M). (a) (the upper left plot) : for the bond activation probability : p cr = 
0.2482 ; (b) (the upper right plot) : for p = 0.245 ~ p cr ; (c) (the lower left plot) : for 
p = 0.35 > per ; (d) (the lower right plot) : fluctuations of the quantity : m = M\ — S m ax , 
where Mi = is the first moment of the fragment-size probability distribution, 

at the percolation threshold p cr are plotted in the scaling variables of delta-scaling (^ with 
(5 ~ 0.80 . For more details, see the description in the text. 
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rapidly, as can be seen in Figs, la and lb (the upper right plot). The latter 
plot is obtained for the bond activation probablity p — 0.245 , very close to 
p cr . Results in Fig. lb are plotted in the scaling variables (j|) for S = 1 , 
even though slight deviations between different calculations can be seen in the 
shoulder region and in the tail for positive zm . Fig. lc (the lower left plot) 
shows the order parameter distribution in the 'liquid phase' (p = 0.35 > p cr ) . 
One finds the second scaling ( |l2|) , in agreement with the analytical derivation. 
Finally, Fig. Id (the lower right plot) shows what happens at the percolation 
threshold (p = p cr ) if instead of P[Nn] , one plots the probability distribution 
of m = Mi — S max , where Mi = kn(k) is the first moment of the fragment- 
size probability distribution, m (= Nf] ), is related in a non-trivial way to 
the order parameter and, in particular, it conserves the singularity of S max . 
P[m] obeys the delta-scaling (||) for 6 ~ 0.80 , what should be compared with : 
1/(t— 1) = 0.84 , given by the analytical argument leading to the delta-scaling 
(^J) . This signature of phase transition disappears for p ^ p cr , i.e., 8 becomes 
equal 1/2. 

As a second example, let us consider the_FIB model uu which exhibits 
the second order, shattering phase transition E3. In this off-equilibrium case, 
analytical derivation of @ and (||), in principle, does not apply. One deals 
in FIB model with clusters characterized by a conserved quantity, called the 
cluster mass. The anscestor cluster of mass N is fragmenting via an ordered 
and irreversible sequence of steps until either the cutoff-scale for monomers is 
reached or all clusters are inactive. 

Each step in this cascade is either a binary fragmentation of an active 
cluster (fc) — > (j) + (k — j) with a fragmentation rate ~ Fj,k-j (a mass 
partition probability), or its inactivation (fc) — * (k)* with an inactivation rate 
~ Ik ■ The order parameter is here the reduced cluster multiplicity Mq/N or 
the reduced monomer multiplicity, both of them closely interrelated. The total 
fragmentation probability pp at each step of FIB cascade is : 

fc-i fc-i 
PF {k) = Fj, k -i (h + £ F^)- 1 . (13) 
i=i i=i 

The cluster mass independence of PF(k) at any step until the cutoff-scale 
characterizes the critical transition region. FIB process is self-similar in this 
regime. For pp > 1/2 , the anomalous dimension (^) (now iV|?7| = Mo) varies 
from to 1, what is different from the limits on in the equilibrium systems. 
The average multiplicity of inacitve clusters is EH : < Mq >^ N T ~ 1 (1 < 
r < 2) , and g = t — 1 . In the shattered phase, the average multiplicity is : 
< M > ~ N , the cluster-size distribution is a power-law with r > 2 and 
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Figure 2: The order parameter fluctuations in the FIB model calculated for two initial sizes 
: N = 2 10 and N = 2 14 , and plotted in the scaling variables of delta-scaling (hi), (a) (the 
upper left plot) : the critical FIB process for pp = 0.875 , a = corresponding to the 
anomalous dimension g = 0.75 ; (b) (the upper right plot) : fluctuations of the quantity : 
m = N — Mo are plotted in the scaling variables of delta-scaling (H) with 6 ~ 0.73 , (c) (the 
lower left plot) : the critical FIB process for pp = 0.7 , a = corresponding to the anomalous 
dimension g = 0.4 ; (d) (the lower right plot) : shattered phase for a = 0, b = —1, r = 4 . 
For more details, see the description in the text. 
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the anomalous dimension is g = 1 . In this phase, pp is an increasing function 
of cluster mass k and, hence, the FIB cascade is not self-similar. Most of 
the interesting physical applications correspond to homogeneous fragmentation 
functions! : F^ j X (N~j) = >? a Fj,N-j , e.g., Fj, N -j ~ \j(N - j)] a . For the 
homogeneous inactivation rate- function : 1^ ~ k b , the critical transition region 
in FIB model corresponds to the line : b = 2a + 1 , and the shattered phase 
corresponds to : b < 2a + 1 . Such homogeneous rates : Fj^-j and J/. , will 
be used in examples shown in Fig. 2. Fig. 2a (the upper left plot) shows the 
scaling function of critical FIB process for pp = 0.875 , a = , what 

yields g = 0.75 . The asymmetry of <l>(z(i)) about Z{\) = , is rather common 
in the critical sector of FIB model til. This sectot, and its characteristic first 
scaling (||) extends to the domain 0<g<l/2tJ, excluded in equilibrium 
models. In this domain (g < 1/2) , the most probable value ofp^(z(i)) is at 
Z(x) = —1 , whereas for g > 1/2 it takes a value close to zn\ = 0O . This can 
also be seen in Fig. 2c (the lower left plot) which exhibits the scaling function 
^(z^) of critical FIB process for pp = 0.7 , a — , for which g = 0.4 . What 
happens if instead of P[M ] = P[Nrf\ , one plots P[N - M ] = P[Nfj] , as 
shown in Fig. 2b (the upper right plot). Analogously as in percolation (see 
Fig. lb), P[Nfj] is scaling (JsJ) with the non-trivial exponent S ~ 0.73 which is 
close to the value 5 = g = 0.75 , obtained using analytical arguments (||). The 
order parameter distribution in the shattered phase (a = 0, b = — 1, r = 4) , 
is shown in Fig. 2d (the lower right plot) in the variables of second scaling 
( |l2| ) for S = 1/2 . Again, one finds an analogy to the situation in the 'liquid' 
phase of percolation (Fig. lc). 



4 Outlook 

The off-equilibrium FIB model shows essentially the same relation between 
criticality and scaling of order parameter fluctuations as it has been derived 
analytically for the second-order equilibrium transitions (Eqs. (||), (||), (|l2|)). 
This is related to the underlying self-similarity which is common to both equi- 
librium and off-equilibrium realizations of the second-order transition. In that 
sense, the scaling laws (Eqs. @, ([l2|)) are the salient features of any 
system exhibiting the second-order transition and the function ^(z^) is a 
fingerprint of the system and its transition. It is therefore important to study 
implications of scaling (||) on statistical properties of the system. Let us first 
postulate the definition of the pseudo-free energy : 



T=-P T 1 \n{< \t)\>P[t)]) 
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(14) 



with a coefficient Pt which is independent of n and characterizes the homo- 
geneous system. If we suppose that the first scaling holds and employing 
the asymptotical definition of the anomalous dimension (Q) , one finds : 

T/N ~ 77 1/(1 - 9 V(^ 1_ff ) , (15) 

what formally corresponds to (^) with (1 — g) instead of both /3/(2 — a) and 
f3/(vS). T appears to be a pseudo-free energy for constrained r\ , and fixed 
both N and fix- This derivation does not suppose that the system is close to 
the second-order transition in thermodynamical equilibrium. In particular, (^|) 
and (0) are expected to hold for any critical behaviour : at the thermodynamic 
equilibrium as we have shown rigorously, but also for the non-thermodynamic 
equilibrium like in percolation, for the ofLcquilibrium final state like in FIB 
model, or for the self-organized criticalitylij. The basic quantity is T . which 
contains thermodynamically relevant information about the analogue of the 
inverse of temperature. The scaling law of T ( |l5|) gives in turn information 
about the critical behavior of the system. Finally, let us remind that the 
formulae (|l^) contains only one exponent g. To have an access to other critical 
exponent, we need to vary the field conjugate to the order parameter. 

To summarize, we have found the new approach to study critical phenom- 
ena both in equilibrium and off-equilibrium systems. This approach is based 
on the existence of universal scaling laws in the probability distribution of both 
the order parameter and its complement, in the second order phase transitions. 
The precise relation between the scaling functions ^(z^) , the nature of or- 
der parameter and the critical exponents yields a new tool for determining the 
combinations of critical exponents even in small systems and for learning about 
the nature of critical phenomenon. We hope, this approach will be useful in 
many phenomenological applications in the strong interaction physics and in 
the condensed matter physics. 
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